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Abstract — In this paper, we derive the exact Cramer-Rao 
bounds (CRBs) for semi-blind channel estimation in amplify- 
and-forward two-way relay networks employing square QAM 
modulation. The derived bounds are used to show that the semi- 
blind approach, which exploits both the transmitted pilots and 
transmitted data symbols, can provide substantial improvements 
in estimation accuracy over the training-based approach which 
only uses pilot symbols to estimate the channel parameters. We 
also derive the more tractable modified CRB which accurately 
approximates the exact CRB at high SNR for low modulation 
orders. 

Index Terms: Amplify and Forward, Cramer-Rao bound, 
Semi-blind Channel Estimation, Square QAM Modulation, 
Two-way Relays. 



I. Introduction 

Amplify-and-forward (AF) two-way relay networks 
(TWRNs) [ 1 1 have recently received a lot of attention as a 
spectrally efficient approach for bidirectional communication 
between two terminals. In previous works on TWRNs, it 
is often assumed that the channel parameters are perfectly 
known at the terminals. In reality, however, these parameters 
have to be estimated, and an estimation error is always 
incurred. Since part of the channel information is required for 
self-interference cancellation at each terminal, the presence 
of estimation error means that self-interference cannot be 
completely cancelled, and the residual interference reduces 
the overall performance of the system. 

The problem of estimating the channel in AF TWRNs has 
been addressed in a number of recent works (2)-||9). Most 
of these works adopt the pilot-based approach to channel 
estimation (c.f. (2j), which imposes an additional burden on 
the system and reduces the overall spectral efficiency. It is 
therefore important to investigate the application of semi-blind 
channel estimation which mitigates this burden by exploiting 
both pilots and the transmitted data to estimate the channel, 
or blind channel estimation which completely avoids this 
burden by relying only on the transmitted data for channel 
estimation. Researchers have recently started investigating 
such approaches (7), (SJ in order to obtain better tradeoffs 
between estimation accuracy and spectral efficiency than the 
pilot-based approach. In (7), a semi-blind approach which 
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jointly estimates the channel and detects the transmitted data 
for MIMO-OFDM TWRNs was developed using the expec- 
tation conditional maximization (ECM) algorithm. In (8|, a 
blind channel estimation algorithm was proposed for TWRNs 
employing constant modulus (CM) signalling in the form of 
MPSK modulation. 

In order to evaluate the potential of semi-blind and blind ap- 
proaches in the AF TWRN context and to provide a benchmark 
on the performance of such algorithms, it is very useful to 
know the Cramer-Rao bound (CRB) on achievable estimation 
accuracy. This bound is more easily obtained for the pilot- 
based approach, as has been done in (2j. However, for semi- 
blind and blind approaches, the task of deriving the CRB 
is more challenging because of the complicated form of the 
likelihood function when the statistics of the data symbols 
are taken into account. In (8j and ||9), the derivation of the 
blind CRB for the case of CM signalling was simplified by 
ignoring the statistics of the data symbols and instead treating 
them as deterministic unknowns. To the best of our knowledge, 
however, the exact CRBs that take into account the statistics 
of the data have not been derived before for AF TWRNs. 
In this work, we fill this gap and we derive the exact CRBs 
for semi-blind channel estimation in AF TWRNs employing 
square QAM modulation. Due to its bandwidth efficiency, this 
class of modulation schemes is very important for today's high 
datarate applications. The derived bounds also cover the blind 
and the pilot-based scenarios as special cases. We use the 
derived bounds to explore the potential of the semi-blind and 
blind approaches and show that substantial performance gains 
over the pilot-based approach can be achieved. 

The rest of the paper is organized as follows. In Section |H| 
we present the system model. Our derivation of the CRB is 



presented in Section III Simulation results are presented in IV 
Finally, our conclusions are in Section [V] 

II. System Model 

We consider the half-duplex TWRN with two source nodes, 
71 and 72, and a single relaying node 1Z, shown in Fig. [T] 
The network operates in quasi-static flat-fading channel con- 
ditions. Each transmission period is divided into two phases. 
In the first phase, 71 and T2 simultaneously transmit to 1Z, 
and in the second phase 1Z broadcasts an amplified version 
of the received signal to both terminals. The semi-blind 
approach employs both pilot symbols and data symbols to 
estimate the channel parameters. More specifically, prior to 
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used in the presence of random nuisance parameters and is 
more tractable than the exact CRB. In deriving these bounds, 
we also take into account the parameter \h 2 \ 2 which, though 
not required for detection, appears in the likelihood function 
and thus affects the estimation performance. The unknowns 
parameters are thus a, b and \h 2 \ 2 and are collected into the 
real vector 9 4 [3?{a}, 3{a}, 5R{6}, 3{6}, \h 2 \ 2 ] T ■ 



Fig. 1. The two-way relay network with two source terminals and one relay 
node. 



transmitting data symbols, each terminal transmits a block 
of L pilot symbols. We denote by ti = \tu, . . . ,£il] T and 
t 2 = [t2i ) • • • > *2i] the vectors containing the pilot symbols 
transmitted by 71 and 72. respectively. The received signal 
vector at 1Z during the training period is r t = + <7i*2 +<*>, 
where hi and gi are the complex coefficients of the flat- 
fading channels 71 — > 1Z and 72 — ^ 72., respectively, and 
u is the circular complex white Gaussian noise with mean 
zero and covariance a 2 I (denoted as CCJV(0, a 2 I)). The relay 
broadcasts Ar t , where A > is the amplification factor. The 
corresponding received signal vector at terminal 71 is 



Ahih 2 ti + Agih 2 t 2 + Ah 2 uj + wi 



(1) 



where u>i is also CCAf(0, a 2 I). 

After transmitting the L pilots, 71 and 72 transmit TV data 
symbols each. We denote by s\ = [su, . . . , siat] t and s 2 = 
[s 2 i, . . . , s 2 n] t the transmitted data symbol vectors of 71 and 
72, respectively. The received signal vector at 7Z is r = hisi + 
gis 2 + n, and the corresponding received signal vector at 71 
is 

z = Ahih 2 Si + Agih 2 s 2 + Ah 2 n + rti (2) 

where n and ri\ are CCJ\f(0, a 2 I). The complex channel 
coefficients hi, h 2 , gi and g 2 are assumed to remain fixed 
during the transmission of the L pilot symbols and N data 
symbols. 

We assume that both terminals employ square QAM mod- 
ulation with possibly different modulation orders and trans- 
mission powers. Without loss of generality, we focus on the 
derivation of the CRB for channel estimation at terminal 71. 
For square QAM modulation, the total number of constellation 
points is M = 2 2p , where p — 1,2,3,.... Denoting by 
d the intersymbol distance and letting d p = |, the set of 
constellation points used by 72 is given by S — {±d p (2i— 1)± 
jd p (2£ - 1)}, i, I = 1, . . . , 2P" 1 |j?0j. The average transmitted 
power at T 2 is P 2 = E {|s 2 fc| 2 } = ^d 2 . We also let P x be 
the avarage transmitted power at 71- Furthermore, we assume 
that the noise variance a 2 is known at 71- In our work, we 
are interested in deriving the CRBs for the estimation of the 
composite channel parameters a = hih 2 and b = gih 2 , which 
are sufficient for data detection. 

III. Cramer-Rao Bounds 

In this section, we derive the CRBs for the estimation of 
a and b. In addition to deriving the exact CRB, we will also 
consider the modified CRB (MCRB) [11] which is commonly 



A. Exact Cramer-Rao Bound 

To derive the exact CRB for semi-blind channel estimation, 
we consider the joint likelihood of z t and z. Let z = 
[zf,z T } T , the joint likelihood is given by 
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The corresponding log-likelihood function is 
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(4) 



where C = a 2 (A 2 \h 2 \ 2 + 1). Let 1(9) be the corresponding 
Fisher information matrix (FIM), and let le^e- be the joint 
Fisher information between the parameters Oi and 9j, where 
i.j = 1, . . . , 5. Furthermore, to simplify our notation we let 
a R — 5R{a}, a-i — 3{a}, bn = $l{b}, b[ = and r = 

\h 2 \ 2 . The matrix 1(9) is given by 



1(9) = -E 
IT 



d 2 C(z-0) 
dOdO T 



a a 

l ab 
tT 



where 



Ibb — 
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h R ,T 



(7) 

Before proceeding to obtain closed-form expressions for 
the elements of 1(9), we will first factorize the likelihood 
function by taking into account the symmetric structure of 
square QAM modulation, following the approach proposed 
in (10). This factorization will make it feasible to derive 
analytical expressions for the elements of 1(9). We begin by 
rewriting the likelihood function in <[3j as 

f(~z-,e) 1 
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(8) 
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Now, we let 

D k (6)= ^ e-c' A2|b|2|s2|2+ ^ 5R{(2:fc " Aaslfc) * f,S2} . (9) 

s 2 €S 

The inherent symmetry of the constellation set allows us to 
write |9]) as a sum over the symbols in the first quadrant. If Qi 
is the set of constellation symbols that lie in the first quadrant, 
S can be partitioned as S = QiU(— Qi)UQ*U(— Q\). Hence, 
we may rewrite D k (8) as: 

D k{0) = e-* A2|b|2|s212 f e ^m^-Aas lk rbs 2 } 
s 2 eQi ^ 

_|_ e ^T$l{(zk-Aas lk ybs* 2 } , g - 2£ 5R{(z fc - Aas lk )* bs 2 } 

(10) 

Noting that 5R{(z fe - Aasi k )*bs 2 } = ^{{z k - 
Aasi k )*b}^R.{s 2 }-'^{(z k -Aas lk )*b}'^{s 2 }, we rewrite ( fTO] ) 
as 

-*A 2 |6| 2 | S2 | 2 



x cosh 
x cosh 



2A 
~C 
2A 
~C 



di{(z k - Aas lk )*b}%i{s 2 } 
- Aas lk )*b}Q{s 2 } 



(11) 



Moreover, for s 2 E Qx, we have that $t{s 2 }, ^s{s 2 } £ {(2i 
l)d p }, i = 1, . . . , 2 P ~ 1 . Hence, (jTTJ becomes 



2 P ~ 1 2 P— 1 

,0,(0) = 4 ^ ^ e -£A°|l,| 2 ((2i-l) a + (2/-l)Xx 
i=l f=l 

"2A 
"2A 



cosh 
cosh 



-(2i-l)dp»{(2 fc -i4<wifc)*6} 
-(2^-l)d p 3{(z fe -^as lfc )*&} 



(12) 



From (12) , we can see that, similar to the case in fTO], D k (6) 
is the product of two terms, one depending on di{(z k — 
Aasi k )*b} and the other depending on — Aasi k )*b}: 



D k (8)=4F e (u k )F e (v k ), 

2 P ~ 1 

Fe (t) 4 J2 e-h^^W^-if CQsh f^P {2t _ 1)t 



(13) 



where 



and 



V c 

u fe = $i{(z k - Aasi k )*b} 

= $t{z k - Aas lk }b R + 3?{z fc - Aasi k }h 

w fe = - Aas lk )*b} 

= ^{z k - Aasikjb! - Q{z k - Aasi k }b R . 



(14) 
(15) 



(16) 



Before proceeding, we simplify our notation by letting /3, = 



A d p 
~G~ 



2i - 1) and 7l = ^S-(2i - l) 2 for i = 1, . . .^P" 1 . 



Using the newly defined 7.;, we can write -Fg(t) as 
2 p-i 

F e (t)^ X! e ~ 7,|b|2 cosh ( 2/3 ^ ■ 

Using (jT3j, the likelihood function becomes 



(17) 
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As we shall see shortly, the RVs u k and v k are independent, 
a fact which will simplify our derivation of the FIM. It is 
also useful to define two new RVs, x k — $t{z k — Aasi k } and 
y k = %{z k - Aas lk }. The pairs {u k ,v k } and pair {x k ,y k } 
are related through the following linear transformation 

(19) 

Both pairs of RVs will be used in deriving the elements of the 
1(0). It is easy to see that the joint PDF of x k and y k is 
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fx.Y{x k ,y k )= e a Fg(b R x k +biyk)Fe(biXk-bRyk)- 

* MC (20) 
Using (jT9]» and ( pO) , we obtain the joint PDF of u k and v k : 
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e ~cW F e (u k )Fg(v k ). (21) 



ttMC|&| 2 

It is clear from ( |2T| that the RVs u& and are independent 
and identically distributed (i.i.d.) with respective PDFs 
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and 



fu{u k ) 
fv(v k ) 
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(22) 



(23) 



a/M^C^ 2 

Going back to the log-likelihood function in we may 
now rewrite it as 

1 

C' 
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C(z;0) 



(N + L) log(^C) - ^||z p - Aati - A6t 2 f - 



N 



N 



log («*) + X! lo S F » ( Vfe ) ' 
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fe=i 



(24) 



From ( |24) we see that the main task in obtaining analytical 
expressions for the elements of 1(6) is the evaluation of the 
expectations E { "^A™ } and E { } for i, j = 

1, . . . , 5. Letting 



B 



d 2 Fg(u k ) 



F (u k ) 



G 



dFgjuk) dF e (u k ) 

{if) A "i d9 l 
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H 



d 2 F e {v k ) 
(ij) a g^j 



ir. 



(ij) a ggj ag 3 



(25) 



(26) 
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we have that 



and 



2 log F e (u k ) _ m _ 



o 2 log F g (v k ) =H (ij) _ w m_ 



09,09, 



(27) 



(28) 



Despite the factorization of the log-likelihood function, 
the derivation of analytical expressions for the elements of 
1(9) requires tedious calculations. Due to space limitations, 
we will provide detailed derivations for only some of these 
elements. For the remaining elements we will provide only 
the resulting analytical expressions. We begin with the first 
diagonal element of the FIM, I aR ,a R - We have 

2 £(z;8) 



E 



N 



2A 2 „ 2A 2 „ 

N 



(29) 
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We show in Appendix [A] that 

E{ B W} = ^=J2PMsUb} 2 , OO) 
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To obtain E •^G^. 11 ' 1 we need the first derivative of Fg(u k ) 
with respect to a R , which is given by 



dFe{u k ) 
0a R 



-2AU{s* lk b} ^ fte-^l fc l 2 sinh[2fau k ], (32) 



i=i 



From ( f32] l, we see that we can evaluate E |g[ U ' ) | using the 
PDF fu(u k ) in d22). We thus obtain 



where 



and 



E{G{ 11) } = n{si k b} 2 r 

00 f(t) 



(33) 



r 8A 2 

1 ^7rMC\b\ 2 J~oc F e {t) 



e °w*dt, (34) 
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/(f) = ^ fo-vM' smbpfiit]. 
i=l 

Moreover, it can be easily verified that 

E{^ n) } = 3K fc &} 2 r l5 



(35) 



(36) 



Hence, 
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+ H. 

O.R — (J 1 (J 
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(37) 



Regarding the second diagonal element of 1(9), it can be 
easily shown that l auai = I aR , aR - 

For the third diagonal element of 1(9), we have 
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Furthermore, we show in Appendix [B] that 

2 p-i 2 p-i 

J2 E mift- 



Ejsi 33 ^ = 2 
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(38) 
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To find E |g^ 33 ' ) |, we obtain the derivative of Fg(u k ) with 
respect to b R , which is given by 



qi{xk,Vk) 



0F g (b R x k + biy k ) 
0b R 



2 p-i 



= - E 2 ^i b R e 7,|b| cosh[2 Pi(b R x k + b I y k )} ^ 

i=l 
2 p-i 

+ 2 2fcx k e-^ 2 wabppiibnXH + b lVk )]. 



Let r 2 4 E {g^ 33 ' } = E { ^ffiffi' }■ Clearly, T 2 cannot be 
evaluated using the PDF fu(u), and we need to use the joint 
PDF fxy(x,y) in ( pO) , which means that double integration 
is required. Using (p0|>, we obtain 



ql{x,y)F e (b I x ~b R y) _^+j L 



ttMCJJ F 9 (b R x + b iy ) 



e c dxdy. 



(41) 



Moreover, it can be easily verified that E | Vt^ 33 ^ j 
E |Gf 3) |, which implies that 



l b R ,bp, 



op-1 op-l 

327V V- V- ,2 
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(42) 



2AT? 



A very similar approach can be followed to evaluate 
thus obtaining 



Ibj,b, 



where 



2A 2 
~~C 



327V 



ft-TEE^ + 2JVr 3 , (43) 
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2 p-i 
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For the fifth diagonal element of 1(9), I T , T , we have 
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It can be easily shown that E {\\z v — Aat\ — ^46*2 1 1 2 } = LC 
and E{||z- Aasi|| 2 } = NA 2 \b\ 2 P 2 + NC. Moreover, we 
show in Appendix [C] that 



E{5f)}=E{ 



i—l ^ 

(47) 

To obtain e|g[, 55 ' ) |, we first take the derivative Fg(u k ) 
with respect to r: 

2P- 1 

q 3 (u k ) 4 ^^1^1 = £ AV 2 A 2 |6| 2 e-^l 6 l 2 cosh[2ft Ufc ] 
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Letting T 4 = E [g^ 55) |, we thus get 
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Clearly, we also have that EjW^ 5 } = T 4 . Using the above 
results, we get 



I T ,r = -(N + L) 
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We now consider the off-diagonal elements of 1(0). First, 
we show in Appendix [D] that / aR ,aj = 0. The remaining 
elements of 1(9) can be obtained using similar approaches to 
those employed so far, and we will only provide the resulting 
analytical expressions. Going back to the sub-block matrices 



of 1(9) in (BJ, the matrices I aa and J&& are given in ( BTj ) 
and (|52"1>, respectively, at the top of the next page, where 
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For the remaining submatrices, we have that 
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Having derived the FIM matrix 1(9), the exact CRBs on a 
and b can be obtained by taking the inverse of 1(9). We havaH 



CRB a = [I(9)~ 1 } 11 + [I(9)- 1 } 22 , 



and 



Ci?B h = [I(6)-%3 + [I(9)- v 
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(60) 



The bounds in ( |59] l and ( |60| l provide convenient benchmarks 
for the performance of estimators of a and b. For N = 0, 
they cover the case of fully pilot-based estimation and for 
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L = they cover the case of blind estimation^] In Section 
we will use these bounds to compare the the semi-blind ap- 
proach with the pilot-based and blind approaches. The obvious 
shortcoming of the exact CRBs is their high complexity. We 
next consider a more tractable alternative, the modified CRB 
(MCRB). 

'The notation is used to refer to the (i,j')th element of the matrix 

A. 

2 We note that in the absence of pilots the channel parameter a remains 
identifiable thanks to the presence of known self-interference symbols. The 
parameter b, however, suffers from an inherent ambiguity and is only locally 
identifiable. The CRB is still defined in this case, as pointed out in |12| . In 
practice, however, pilot symbols are required to resolve the ambiguity. 
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B. The Modified Cramer-Rao Bound 

Due to its tractability, the MCRB is commonly used in the 
presence of random nuisance parameters fTT) . In our case, the 
nuisance parameters are the data symbols s 2 k, k = 1,...,N. 
To obtain this bound, the data symbols are initially treated 
as deterministic unknowns, which greatly simplifies the like- 
lihood function and the resulting FIM. The statistics of the 
data symbols are then taken into account by averaging the 
FIM over the data symbols. We call the resulting matrix the 
modified FIM (MFIM) and denote it by J. It can be shown 
that the matrix J is given by ( |oT| on top of the next page. 
Denoting by MCRB a and MCRB the resulting bounds for 
parameters a and b, respectively, it can be shown that 



MCRB a = 

and 

MCRB b 
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The above expressions are much simpler and more tractable 
than those in |59| and ( |60] >. However, as we shall see in 
Section |IV] they are loose and do not reflect the impact of 
the modulation order on the estimation accuracy. However, 
for low modulation orders they are close to the true bounds at 
high SNR. 

IV. Simulation Results 

In this section, we use MATLAB simulations to investigate 
the behavior of the derived CRBs for a and b for the cases 
of blind, semi-blind and pilot-based estimation. All plots are 
generated using the Monte-Carlo approach and are averaged 
over a set of 100 independent realizations of the channel 
parameters hi, h 2 , g\ and These realizations are generated 
by modelling hi and h 2 as correlated complex Gaussian 
random variables with mean zero, variance 1, and a correlation 
coefficient g = 0.3. Similarly, we model gi and g 2 as 
correlated complex Gaussian random variables with the same 
mean, variance and correlation coefficient, but independent of 
hi and h 2 . To generate correlated complex Gaussian random 



variables we follow the approach proposed in [13 J. In order 
to see the effect of the modulation order on the CRB, we 
consider four modulation orders, M = 4, M = 16, M = 64 
and M = 256 in all our plots. In all our simulations, the pilots 
are generated using M = 4, and the pilot vectors of the two 
terminals are are orthogonal to each other. The SNR is defined 
as 10 log =9f. 

We begin by comparing the CRBs of the semi-blind ap- 
proach and the pilot-based approach. In Figs. [2] and [3] we 
plot versus SNR the semi-blind CRB and the pilot-based CRB 
for parameters a and b, respectively. The number of pilots 
is L = 8, and the number of transmitted data symbols is 
N = 32. We also plot the corresponding MCRB for M = 4 
and M = 256 in both figures. As we can see from both figures, 
the semi-blind CRB is substantially lower than the pilot- 
based bound for all 4 modulation orders. This shows that, by 
making use of the transmitted data symbols in addition to the 
pilots, the semi-blind approach can provide substantial gains 
in estimation accuracy over the pilot-based approach, thus 
providing a superior tradeoff between accuracy and spectral 
efficiency. We also see that the lower the modulation order 
the higher the achievable accuracy, with the best accuracy at 
M = 4. In addition, Figs. [2] and [3] show that the MCRB is 
generally loose compared to the true CRB and is not sensitive 
to the modulation orders, and thus does not reflect the effect 
of the modulation order on the estimation accuracy. However, 
MCRB provides a good approximation of the true CRB at 
high SNR for M = 4 and M = 16. 

In Figs. [4] and [5] we plot the semi-blind CRBs of a and 
b versus N, respectively. The number of pilots is L = 8, 
and the pilot-based CRB is plotted as a reference. As we 
see from both plots, as N increases, the semi-blind approach 
can provide increasingly better accuracy compared to the 
pilot-based approach. The sample size N is constrained by 
the coherence time of the channel during which the channel 
parameters a, b remain fixed. Hence, the longer the channel 
coherence time the more attractive the semi-blind approach 
becomes. 

We next compare the CRBs of the semi-blind approach 
and the blind approach. In Figs. [6] and [7] we plot versus 
SNR the semi-blind and blind CRBs for parameters a and b, 
respectively. For the semi-blind case, L = 8 pilots and = 32 
data symbols are employed. For the blind case, N = 40 data 
symbols are used. We see from Fig. [6] that the semi-blind 
approach provides better accuracy for the estimation of a. As 
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the modulation order increases, the gap in favor of the semi- and 

blind approach becomes more significant. In Fig. [7] however, , 2 2P ~ 1 2 

we see that the blind CRB for b completely deteriorates at E \ B k } = ^7^^{ s *k^} 2 E ^ 69 - ) 

low SNR, reflecting the difficulty in the estimation of b in the i=1 

absence of pilots. The difference between the behaviors of the Using very similar steps, we can also show that 

blind CRB for a in Fig. [6] and for b in Fig. [7J reflects the fact 

that in the absence of pilots the estimation of a is an easier mo 8A 2 2 2 

task than the estimation of b due to the presence of the known ft i J = 

self-interference symbols. ' _J 



V. CONCLUSIONS APPENDIX B 



In this paper, we derived the exact CRBs for semi-blind In this appendix, we prove ([39Jl. The second derivative of 

channel estimation in AF TWRNs employing square QAM F g (u k ) with respect to b R is given by 

modulation. As an more tractable alternative, we also de- x 

rived the modified CRB. Using Monte-Carlo simulations, we d 2 Fg(u k ) _ -sr-^ . 2 ^ 2 2 ) - 

showed that the semi -blind approach can provided substantial db% ^— ^ R 



•7«i&r 



gains over the pilot-based approach. The semi-blind approach 
also seems to be more attractive than the blind approach, 



x cosh[2Pi(b R x k + 6/2/fe)] 



which deteriorates at low SNR. For future work, it remains to ^ (71) 

design efficient low-complexity semi-blind algorithms whose _ 2-* & ^i^R x k e 7 ° smh[2 (3 i(b R x k + 

performance approaches the semi-blind CRB. 



i=l 



Appendix A 

In this appendix, we prove ( pO) and ( |3T| . Recalling that 

u k = 3?{^fc — Aas lk \b R + — Aasik}bi, the second 

derivative of Fg(u k ) with respect to or is given by 



2 4f3fx 2 k e-^ 2 coah[2^{b a x k + b lVk )]. 



i=i 

We let 

2 p-i 



= 4A 2 K{^ fc &} 2 £ ft 2 e-^l 2 cosh[2ftu fc ]. 



i=l 
2 p-i 

(64) = - ^ 8/3 i 7 l 6i ? a; fc e" 7i|61 sinh[2A(&fl^fc + &ry*)] 

i=l 

2 P-1 

f iM9fl . n r 3 -E 4 ^ e ~ 7l|fc|2cosh [ 2 ^(^ + ^)]- 



1=1 

Thus, 2 p_1 



I ^K) J ' (72) 

(65) We next indEj^), eI-et^), and E ( -gr^r ) . Us- 
ing d66|, we obtain 



2=1 

Using the PDF for in ( |22| i, we obtain 

E < L /* fc H = -— /cosb-2ftie ^dt. T t 1 , , , 

I F e (u k ) / ^fMCW^ E =^TtE( 8 ^- 4 ^)- (73) 

(66) 

Moreover, for a > 0, it can be verified that We next consider E j F ^ tfc) j, We have 



V— dt = v /|c4. (67) _ f T2 ^ 



Hence. lFeK)/ S (74) 

f cosh[2^iU fe ] 1 = 2 | 6 |2 f x fc sinh[2/3. t (b J? x fc + b f y fc )] 1 

~' *K«fc) J \/m I Fg^xV+biyk) J 
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Using the PDF fxx{x,y) in {20) , we obtain 



E 



-X 



x k smh[2(3i(b R x k + biy k )] 

Fg(b R x k + b lVk ) J itMC 

x smh{2(3i(bB.x + bjy)]Fg(bjx — b R y)e o c?xdy 



After some calculations (the details are skipped for brevity) 
we obtain 

x l cosh Wi ( b R x k + b lVk) 



E 



ttMC 



Hence, 



E 



= 7.|f»| 2 



r 3 

F e {u k ) 



Fg{b R x k + biy k ) 
C 2C 



M VM 



> 4C — > . > 



(81) 



IP — 1 



E^ + tt^ E^ 2&; 



9 p— i 



-E 



-7*161 



«=1 



-oo 

oo 

+ / Xe~ 3 £ + (2Pit>R-2&lbi)x dx j e -£ + (2p l b I +2/3 i b R )y dy 

— oo 
oo 

:z . e -#-(2/3 I fcfi-2ftfc / )x d;z , J e -£—(2Pibi+2Ptb R )y d y 

— oo 

OO 

xe -^-(2Pib R +2/3 e bj)x dx f e -^--(2p i b 1 -20 e b R )y dy 



16 

M 



i=l 
2 p-i 2 p-i 



i=l 
2 



(82) 



E E 7*7^{6} 5 

i=l f=l 

Adding {73), {78) and {82), we finally obtain 

2 p-i 2 p-i 



i=l £=1 



(83) 



Following very similar steps, it can be shown that 

E{fl"i 33) }=E{B. (33) V 



(75) 

In the above expression, the single integrals in x can be eval- 
uated using the following result which holds for a > [14]: 



Appendix C 

In this appendix, we prove {47). Taking the derivative of 
Fg(u k ) twice with respect to r, we obtain: 

d 2 F e {u k ) ~ ' 
<9t 2 



te 



-af*-26t 



dt = —\ I —Se <* 



(76) 



Y^Ua^f\b\ 4 -^f3f\b\ 2 ) 

x e - 7l|b|2 cosh[2/3 4 u fc ] 

E ( 4 4^A - 4^#|6| 2> ) ^^ 7! "" 2 Binh[2A« fc ] 



while the single integrals in y can be evaluated using {67). 
After evaluating all the single integrals in {75) , we finally 
obtain 

' xsinh[2^(b R x + biy)}\ _ 4C op _ 1 



4 E 4=J^ v fc e " 7 *'" 2 cosh [ 2 ^^] 



(84) 



E 



F e (b R x + biy) 



M 



2 P-l e 7i|6| (77) Since B (55) = 



Therefore 



to evaluate the terms E 



it is clear from <f84b that we need 



E 



T 2 



Fg{u k ) 

We next consider E \ „ T 3 , > 

op- 1 

T 3 



16 



E ^ 



(78) 



f cosh[2ft« fc ] "I w f M fc sinh[2fttt fc ] \ , 
{ F e (u k ) Fe (u k ) / and 

{^H- The first term is found in ^ For the 



1=1 
We have 



E 

second term, we have 

u fc sinh[2/3iU fe ] 



E 



9{u k ) 



y/-KMC\b\< 



E 



Fg(u k ) 



E 4 # 



2 p -li\bV 



isinh[2/3 4 i]e d^i 2 c/t. 



(85) 



E 



i=l 
^2 



(79) 



gcosh[2 / 8 i (b fl Xfc + b iyk )] \ 
■ J 



The above integral can be evaluated using the result in {76 
thus obtaining 



Fg(b R x k + biy k ) 
To evaluate {79), we can follow the same approach that we 
used in {75) . However, instead of using {76) , we would use 
the following result, which holds for a > |14|: 



E 



u k smh[2/3jU k } 

F 0(u k ) 



'M 



t 2 e~ at2 - 25t dt 



7T r2 5 2 1 f ( 2 



For the third term, we have 

_ j u| cosh[2/3 l w fe ] 1 2 



(80) 



C|6| 2 e 7 ' |b|2 A- (86) 



t 2 sinh[2ftt]e ci"i 2 di. 

(87) 
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We evaluate the above integral using ( |80] >, obtaining 

1 



We next consider e{g[ 12) } and e{V^ 12) }. To obtain 



._ j u\ cosh[2ftttfc] 



Fg{u k ) J VM 

Using ([68]), ([86]) and ([88]), we finally see that 

2 p-i 



e^ b \ (2C 2 \b\ 4 /3 2 + C\b\ 2 ). 

(88) 



e|g^, 12 ' ) | we need 9F Q^ k ^ which is given in ([32]) and we 

dF e (u k ) 



also need — ^ k> , which is given by 



E 



K 5) } = E 



2A A a 4 



(89) 



dF e (u k ) 
dai 

Hence, 



2 p-i 



= -2AZ{s* lk b} fte^'l'l sinh[2Au fc ]. (95) 



G>F„(u fc ) a-F e (n fc ) 



The above derivation can be replicated to obtain E{i?| 5 } 
' ~'~~> u k with Vk, which means that E |i?^ 55 H = 



/^i(12) dan dai 

k ~ F e (u k )* 



e{b. (55) \ 



2 P — 1 

/ J2 ftsinh[2ft Ufc K 2 
4A 2 3?{s* fe 6}3{si fc 6} ( — I 



Appendix D 

In this appendix, we show that E j g a X/ [ = 0. We have 



(96) 



For W k l2 \ we need dF g^" k ^ and dF g^" k \ which are respec- 
tively given by 



E 



(d 2 £(q;6)\ ^ f ™ n 2 ) \ a „ , , 2"- 1 

(^r}=E E H } } ^ = -2A^ fe& }Efte-l^ S inh[2 W! (97) 

*=i (90) 9a fl 



JV 



and 



fc=i 



To prove that E j d d ^^^ j = 0, we will show that 

E{Bf )} = -E{flP>} and E {<#»>} = -e{^}. 
We have 



dF e {v k ) 



2 p-i 



9a/ 



2ASi{sJ fe &} fte" 7 * 161 " wahp&vk]. (98) 



Thus, 



3*8 K) dF g (v k ) 



2 P-1 



(12) dan dai 



9 -^0 = 4A 2 n{ S l k bm S l k b} £ ft 2 e 



-7 s |b| 2 



i=l 

x cosh[2/3i?/fe] 



(91) 



2 p- 



and 



d 2 F e (v k ) 



da Rdai = -"'WUbVH'hfi} g A 2 e-- Ib ' 2 (92) 

x cosh[2/3jffe] 



Hence, 



E 



{s[, 12) } =4A 2 5R{ s * u ,6}3K fc 6} ]T A 



2 e -7.ifcr 



x E 



COsh[2ftM fc ] 



(93) 



and 

E 



{^ 12) } = - 4A 2 Vl{ S * lk b}Z{ S * lk b} E A 2 e-^ |b|2 

i=i (94) 

rcosh[2^]| 

Since Ufc and i!& are i.i.d., it is clear that IE ^ i?£ 12 ' 1 j> = 

-eK 2) }. 



/ £ ft sinh[2ft Vfc ] v 2 

= -^i^m^fi}[^^ — J • 

(99) 

Since u k and v k are i.i.d., it is clear from ([96]) and ([99]) that 
E { G i 12) } = - E {^fc 12 '}' which completes the proof. 
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Fig. 3. Semi-blind and pilot-based CRBs for the estimation of b plotted 
versus SNR for Af = 32 and L = 8, and for M = 4, 16, 64, 256. We also 
plot MCRB b for M = 4 and M = 256. 



Fig. 2. Semi-blind and pilot-based CRBs for the estimation of a plotted 
versus SNR for N = 32 and L = 8, and for M = 4, 16, 64, 256. We also 
plot MCRBa for M = 4 and M = 256. 
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Fig. 4. Semi-blind CRB for the estimation of a plotted versus N for L = 8, 
and for M = 4, 16, 64, 256. The pilot-based CRB is shown as a reference. 
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Fig. 5. Semi-blind CRB for the estimation of b plotted versus N for L = 8, 
and for M = 4, 16, 64, 256. The pilot-based CRB is shown as a reference. 



or 
o 



10 



10" 



10° 



10 




10 



- Semi-blind, M=4 

- Semi-blind, M=16 

- Semi-blind, M=64 

- Semi-blind, M=256 
Blind, M=4 

Blind, M=16 
Blind, M=64 
Blind, M=256 




10 



15 
SNR 



20 



25 



30 



Fig. 7. Semi-blind and blind CRBs for the estimation of b plotted versus 
SNR for M = 4, 16, 64, 256. We use N = 32 data symbols and L = 8 
pilots for the semi-blind case, and we use N = 40 data symbols for the blind 
case. 
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Fig. 6. Semi-blind and blind CRBs for the estimation of a plotted versus 
SNR for M = 4, 16, 64, 256. We use N = 32 data symbols and L = 8 
pilots for the semi-blind case, and we use N = 40 data symbols for the blind 
case. 



